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@ Foreign matter detector. 

(g) A foreign matter detector includes an oscillator (1 ) coupled 
to an electromagnetic transducer (10) which includes an 
excitation coil (3) connected to the oscillator (1) and two 
interconnected detection coils (4,5) adapted to produce a 
differential signal therebetween when an object is passed 
through the transducer. A first detector (7A) produces an 
analog signal corresponding to the component of the differen- 
tial signal in phase with the oscillating signal, and a second 
detector (7B) produces a second analog signal corresponding 
to the component of the differential signal in phase with a 
phase-shifted version of the oscillating signal. The analog 
signals are converted to digital signals. As the object passes 
the transducer, the output of the transducer changes in a first 
sense during a first part of the passage and in the opposite 
sense during the remaining part of the passage. A photosensor 
(50) generates a passage signal which lasts for one part of the 
object's passage through the transducer. The passage signal Is 
used to define the digital values of the two digital signals from 
which representative values X and Y are selected. These values 
are used in one or more discriminant equations to determine 
whether foreign matter has been detected. Coefficients for the 
discriminant equations are derived by passing sample objects 
through the detector. The coefficients can be stored, together 
with a code associated with the type of object being detected, 



for later retrieval. 
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Descripti n 

FOREIGN MATTER DETECTOR 



■5 ' Background of the Inv ntion 

This invention relates to a detector for detecting in food products, pharmaceuticai products, etc. the. 
presence of materials such as metais, which, must not be present in the product and which have a different 
quality from the inherent or standard quality of the actual product. 

10 For manufacturers of food products, etc., the presence of foreign matter in raw materials could cause 
damage to processing equipment. Also, the presence of foreign matter in finished products poses problems 
with regard to health and safety of the products. For this reason, detectors have long been used for inspection 
at the point where raw materials for such products are supplied to processing equipment and for Inspection at 
the shipping point of wrapped finished products in order to find any foreign matter. 

15 This type of conventional detector includes an electromagnetic transducer as shown in Fig. 19 or 20, which 
is essentially the same as a differential transformer. The transducer has a primary or excitation coil 3 
generating an alternating electromagnetic field when a high-frequency electrical current is supplied, and a pair 
of secondary or detection coils 4 and 5 electromagneticaliy connected to the coil 3. 
The secondary coils 4and 5 are wound jn opposite directions from each.other and interconnected in series 

20 so that they operate differentially. The voltages or currents induced in the secondary coils cancel each other 
out or offset each other. This forms means for detecting the difference in the signals induced in the secondary 
coils. 

In Fig. 19, the primary coil 3 is wound between and coaxially with the secondary coils 4 and 5. In Fig. 20, the 
primary coil is wound in parallel planes with the secondary coils. 
25 The transducer utilizes the following phenomena which occur when an object being inspected is passed 
through the magnetic field in the direction 31. 

If the object contains iron, the magnetic flux density will increase. The induced voltage of the secondary 
coil 5 adjacent the passing object will first become higher than that of the other secondary coil 4. If the object 
contains a nonferrous metal, the eddy current occurring in the metal will cause a loss of the lines of magnetic 
30 force. As a result, the induced voltage of the secondary coil 5 will first become lower than that of secondary 
coil 4. 

Some food products may contain water or salt which may generate a relatively large signal, for reasons 
similar to those in the case of nonferrous metal, even though the products contain no foreign matter. This 
signal represents the product characteristics of the material, which are called "material effect". 
35 If ferrous or nonferrous foreign matter is contained in these materials, the combination of the signal from 
this foreign matter and the product characteristics result in a composite signal. 

If the foreign matter is small, the signal resulting from it will be small, and there is little difference between 
the signal generated by the product characteristics of the product itself and the combination of this signal and 
that caused by the foreign matter. This makes it difficult to detect the foreign matter. 
40 It is in consideration of problems such as this that detectors such as that described in Japanese Patent 
Provisional Publication SHO.57-198880 have been proposed. 

In the publication, the oscillation signal from an oscillator is supplied via a phase shifter to a transducer 4, 
which outputs a differential signal when an object passes through the transducer. This output signal is 
detected by a first chopper synchronized with the oscillation signal to obtain an in-phase chopped signal, and 
45 by a second chopper having a phase difference of 90 degrees from the oscillation signal to obtain a quadrature 
chopped signal. Each chopped signal Is filtered and independently compared by a level comparator. 

With this detector, to detect foreign matter in products having material effect, by adjusting one of the 
choppers to the phase at which the influence of the material effect is minimized (the phase 90 degrees from 
the phase angle at which the material effect Is maximized), a detection signal with good sensitivity to the 
50 foreign matter can be obtained. 

However, in order to reduce the material effect value to a level one-tenth or less that of the maximum value, it 
is theoretically necessary to keep the allowable error angle of the phase adjustment, within approximately 6 
degrees from the optimum phase angle. In other words, for the actual detection of foreign matter, in order to 
limit the material effect value to a level which does not influence practical use, for example, to a level within 3% 
55 of the maximum material effect value, It is necessary for the allowable error angle of the phase adjustment to 
be within approximately 1 .8 degrees. Thus, precise trial-and-error adjustments are required in order to detect 
foreign matter in products having a large material effect. 

In addition, for the level setting of the level comparator which will be explained later, it is necessary to adjust 
the sensitivity or gain to obtain an appropriate level corresponding to the phase adjustment just described. 
60 Thus, for this detector, the operator is required to perform skilled trial-and-error adjustments for the setting of 
the two variables of the phase adjustment and the sensitivity adjustment. 

In addition, even if the detector is adjusted as described above, it is necessary to give a large reference level 
in advance for offsetting the material effect for foreign matter producing a signal in the same direction as that 
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of the material effect. Consequently, in this phase direction, the detector is used in a condition in which the 
sensitivity has been sharply reduced. Consequently, less detection effectiveness can be expected for metals 
and other foreign matter in this phase direction. 

• ' 5 

Summary f th Inv ntion 

In consideration of the conditions described above, the objective of this invention Is to propose a detector 
for detecting the presence of metals and other foreign matter which does not require the operator to perform 
skilled trial-and-error adjustments such as the two factors of the phase adjustment and the sensitivity 
adjustment, and which is capable of detecting with good-sensitivity foreign matter having a phase direction the 
same as or close to that of the material effect. 

A detector according to this invention comprises: 
means for providing a two-dimensional pair of analog signals representative of electromagnetic parameters of 
an object; . 

means for converting said analog signals respectively into two series of digital values; 
means for determining a pair of representative values each among one of said series of digital values; and 
means for comparing said representative values with reference values. 

20 

Brief Description of the Drawings 
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Fig. 1 is a block diagram of a detector according to the invention; 

Fig. 2 is a block diagram of each of the weighting circuits in the detector; 25 
Figs. 3a-3c show analog signals at various stages in the detector; 
Fig. 4 graphically represents digital values from one of the A/D converters in the detector; 
Figs. 5a-5c graphically represent various weighting functions; 
Fig. 6 graphically represents weighted digital values from one of the weighting circuits; 
Fig. 7a also shows digital values, but obtained when two objects are passed in succession through the 30 
transducer in the detector; 
Fig. 7b also shows weighted values, but obtained by weighting the values of Fig. 7a; 
Fig. 8 is aflowchart of a program for determining primary values; 
Fig. 9 shows monitoring of the weighted values to obtain a primary value ; 

Fig. 10 is a partial schematic view of a modified detector according to the invention; 35 
Fig. 1 1 is a block diagram of a circuitry used with the detector of Fig. 10; 
Figs. 1 2a and 1 2b show signals at different points in the modified detector ; 

Fig. 13 is a two-axis, four-quadrant coordinate system showing vectors which represent the detected 
values and types of metal; 

Fig. 14 is a two-axis, four-quadrant coordinate system representing data for objects having material 40 
effect; 

Fig. 15 is a two-axis, four-quadrant coordinate system representing data for objects having no material 
effect; 

Fig. 16 is a two-axis, four-quadrant coordinate system showing the relationship between the correlation 
coefficient and the region determined by a discriminant equation according to the invention; 45 

Fig. 17 is a two-axis, four-quadrant coordinate cystem showing the region determined by another 
discriminant equation according to the invention; 

Fig. 18 is a flowchart of a program for determining and indicating whether the primary values are in or 
outside a predetermined region; 
Figs. 1 9 and 20 schematically show conventional electromagnetic transducers. 50 
Corresponding reference characters indicate corresponding parts throughout the several views of the 
drawings. 

Detailed Description of Preferred Embodiments 55 

In Fig. 1, an oscillator 1 generates a basic sine-wave electrical alternating current signal having a high 
frequency which is selectable within a range of 30 to 400 kHz. 

An electromagnetic transducer 10 is essentially the same as the conventional transducer shown in Fig. 19. 
and may otherwise be constructed as shown in Fig. 20. SO 

The transducer 10 has a primary coil 3 connected through a first terminal with the output of amplifier 2, and a 
pair of secondary coils 4 and 5 for providing a differentia! signal XY representative of electromagnetic 
parameters of an object being evaluated. 

A photoelectric sensor, which may Include a light emitter 26A and a light receiver 26B, is located in front of 
the secondary coil 5 to detect the passage of the object into the transducer 10. . 65 
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The output line of secondary coils 4 and 5 Is connected through a second amplifier 6 with first inputs of a pair 
of multipiier/detectors 7A and 7B, which are essentially the same as the four-quadrant wave detectors or 
product detectors used in FM receivers. 

The wave detector 7A has a second input connected with the output of first amplifier 2 to receive the 
5 amplified signal as a reference signal in phase with the basic oscillation signal. 

The detector 7B has a second input connected through a phase shifter 8 with the output of first amplifier 2 to 
receive a reference signal out of phase, which is 90 degrees delayed from the basic signal. 

The outputs of detectors 7A and 7B are connected through pairs of smoothers/filters 9A and 9B, 
sample-and-hold circuits 11 A and 11B. analog/digital converters 12A and 12B and weighting circuits 14A and 
10 14B with the input interface 15A of a control or computing device 15. 

The input interface 15A is connected with the light receiver 26B and a device 24 for inputting a region 
coefficient as explained later. 

The control device 15 includes a microprocessing unit (MPU) 15B and an output interface 15C, which is 
connected with an alarm 19, such as a lamp or buzzer emitting a warning indication when foreign matter is 
15 detected. The interface 15C is also connected with an R/q> indicator 20 indicating the magnitude and phase 
angle of the data. 

A timing control circuit 16 is connected with the sample-and-hold circuits 1 1 A and 1 1 B, A/D converters 12A 
and 12B, weighting circuits 14A and 14B and control device 15, to operate them in synchronization. 
A program control circuit 17 is connected with the control device 15, an execution (ON) button 21 which 
20 initiates operation of the control device 15, and a stop (OFF) button 22 which terminates it. 

». .„ The program control circuit 17 is also connected with a selector switch 18 for selecting one of the modes 
"sample test", "equation determination" and "operation" and another switch 23 for selecting one of the sample 
modes "M" and "N" explained later. 
The program control circuit 17 is provided with a keyboard 25 to store data for a type of object and use It for 
25 subsequent inspections. The keyboard 25 has numeric keys for specifying a code number to label the type of 
object after collecting the data, a command key Mi for storing the data in the memory of MPU 15b, and a 
command key Mo for recalling It. 

As shown in Fig. 2 (and the applicant's European Patent Appln. No. 89303673.1 filed on April 13, 1989), the 
weighting circuits 14A and 14B each comprise a predetermined number of stages of delay elements 141, 
30 which are connected in cascade or series, and which may be a shift register. 

The inputs of elements 141 and the output of the last element are connected through coefficient multipliers 
142. respectively, with a summer 143 such as a counter. 

The circuit of Fig. 2 is a non-recursive type circuit, and may be what is commonly known as an FIR (finite 
impulse response) filter. The weighting circuit may be replaced by, for example, a 16-bit microprocessor or a 
35 commercially available multiplier/accumulator designed for digital computations. 

SIGNAL DETECTION AND WEIGHTING 

With reference to Fig. 1, when an object containing iron, for example, is passed for inspection by a belt 
40 conveyor (not shown) through the transducer 10, the electromagnetic induction formed between the primary 
coil 3 and the secondary coils 4 and 5 changes as the object passes. An induced voltage XY in a waveform as 
shown in Fig. 3a is generated as the differential signal from the transducer 10. 

The differential signal XY is amplified by the second amplifier 6, and-then detected by the wave detectors 7A 
and 7B by multiplication by the in-phase and out-of-phase reference signals, respectively. The detectors 7A 
45 and 7B output signals Xde and Yde in waveforms as shown in Fig. 3b. 

The detected signal Xde corresponds to the component of differential signal XY which is in phase with the 
basic oscillation signal, while the signal Yde corresponds to the component of signal XY which has a 90 degree 
phase delay from the basic signal. 

The detected signals Xde and Yde are filtered by the associated smoothers/filters 9A and 9B to be signals Xf 
50 and Yf with waveforms as shown in Fig. 3c, which are temporarily stored in the corresponding 
sample-and-hold circuits 11A and 11B. 

The stored analog signals are converted to digital signals Xdi and Ydi by the corresponding A/D converters 
12A and 12B at regular intervals, and input to the corresponding weighting circuits 14A and 14B, which 
perform weighting processing characterized by a time distribution as will be explained in detail below. 
55 Referring to Fig. 4, the signal Xdi, for example, as a discrete series of digital values x(0), X(T) ... X((n-1)T) is 
input to the weighting circuit 14A at sample time intervals T. (For convenience, this explanation uses zero as 
the time of the first value X(0).) 

Each weighting circuit 14A, 14B has a pulse transfer function for weighting which is a series expressed as a 
polynomial relating to Z" 1 to the order n-1, the number n being odd, as shown in the following Equation (1). 
60 H(z) » Ko + Ki2T' + K2Z-2 + ... Kn.iZ-^1 (1) 

Referring to Fig. 2, the number n of multipliers 142 are assigned weighting coefficients Ko to Kn-1 which are 
distributed in the form of an odd function, as shown in Fig. 5a, 5b or 5c. The coefficients have a 
point-symmetrical distribution, and the central coefficient for the central multiplier 142 is located at the origin 
point of the odd function. 

65 The triangular distribution of Fig. 5a is the reverse form of a typical waveform, which resembles Fig. 4,. 
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representing the signal change with time when iron passes through the transducer 10. The number n is 
determined experimentally according to the duration of the typical wave. The distribution of Fig. 5b takes the 
form of a doublet sine wave, and that of Fig. 5c is rectangular. 

The duration of function H(z) is most preferably substantially the same as the period for which a typical 
material, such as iron, passes through the transducer 1. 5 

When the digital signal Xdi in the form shown in Fig. 4 is input to the weighting circuit 14A, a time series of 
weighted signal values Xw as shown in Fig. 6 are output. Weighted values Yw are obtained in the same way 
from the weighting circuit 14B. 

Thus, a digital signal with the polarities shown in Fig. 4 input to the weighting circuit 14A or 14B will produce 
a weighted signal as shown in Fig. 6, which has the following characteristics: 10 

(1) The amplitude of the maximum value is about twice that of the minimum value; 

(2) The polarity of the minimum value is the opposite of that of the maximum value; this minimum value 
appears twice at approximately an interval t1 preceding and following the maximum value; and 

(3) The data of the output signal which precedes and follows the central maximum value Is substantially 
symmetrical around that maximum value. 15 

Conversely, if the polarities of the input signal Xdi or Ydi are the reverse of those shown in Fig. 4, an output 
signal Xw or Yw having the exact opposite polarities from those in Fig. 6 is obtained. 

In other words, the output signal from the weighting means, to which an input signal such as shown in Fig. 4 
is applied, has the characteristics of a discrete time-series output signal, which is distributed in the form of an 
even function with the origin point having the maximum value, if the time at which the time-series weighted 20 
signal has the maximum value is defined as the origin point of time. 

As shown in Fig. 4, the digital signal Xdi, for example, has a form with two polarities, positive and negative (or 
negative and positive), for a single object * 

As shown in Fig. 7a, if two products containing foreign matter are passed through the transducer 1 one after 
another in close proximity, the signal Xdi will have four peaks P1'-P4'. In this case, a pair of peaks PY and P2' 25 
and a pair of peaks P3' and P4' should be identified separately for the respective products. 

If the signal Xdi were not weighted or otherwise processed specially, the separate signal groupings could 
not be identified. Thus, the peaks P2' and P3' might be incorrectly identified as one pair, due to the polarities 
being reverse. This would result In identification of a different material by means of the computed data and the 
R/o indicator 20. 30, 

As shown in Fig. 7b, the weighted signal Xw can be correctly Identified for the different products, because 
the peaks P1" and P3" correspond respectively to the pair PV-P2' and the pair P3'-P4' in Fig. 7a. Thus, even if 
there is a continuous flow of products, the signals can always be identified for the respective products. 

PRIMARY VALUE DETERMINATION 35 

The MPU 15B has a program, as shown in Fig. 8, to determine or recognize "primary values" out of the 
weighted signals Xw and Yw. The primary values represent the electromagnetic characteristics of an object 
being inspected. 

In order to determine the primary values, It is necessary to monitor a typical change with time of the 
weighted signal Xw, for example, as shown in Fig. 6. The monitoring span should be three to five times the time 40 
interval t1 of the two peaks which occur when foreign matter generates a signal, as shown in Fig. 3c or 4. 

As shown in Fig; 9, an appropriate number 2n + 1 of weighted values Xw obtained during a specified time 
span are constantly monitored. As shown in the following Table 1 , the values Xw are stored in order from the 
newest data to the oldest data as XI to Xi-2n in a series of memory addresses 0-2n. 

45 

Table 1 / - 

Address 0 1 2 ... n ... 2n 

Data Xi Xi-1 Xi-2... Xi-n ... Xi-2n 

50 

When new data is input to address 0, the previous data shift in order from the left to the right, and the oldest 
data at address 2n is discarded. At each time, the data Xi-n at the central address n is monitored to check 
whether or not its absolute value is the maximum, or one of several maxima, among the values Xi to Xi-2n at 
that time (steps 1-2 in Fig. 8) 

If It Is. the following equation (2) is computed to check whether or not there is an identical maximum value 55 
Xi-n + a in the specific vicinity, as generally indicated by the circle A in Fig. 9, of the central address n (step 3). 
O < I a I < 8, Xi-n + a = Xi-n ? (2) 

Provided that the value defines the limit of the vicinity to the central address. 

If there is not, Xi-n is determined as the primary value Xp. If there is, Yi-n and Yi-n+a which have the same 
addresses as those of Xi-n and Xi-n + a are compared. If Yi-n £ Yi-n + a (step 4) , Xi-n is determined to be the 60 
primary value Xp, and the corresponding value Yi-n is determined to be the primary value Yp for the weighted 
signal Yw. 

Conversely, if Yi-n+a is found to be larger than Yi-n, Xl-n is not determined to be the primary value. If this 
occurs, when Yi-n+a is shifted to the central address ri later, that value becomes the maximum value and is 
thus determined to be the primary value Yp. 65 
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If the objects are limited to solid objects or individual batches such as packaged products, because the. 
signals are originally detected Individually, the primary values can be obtained without weighting. 

For example, as shown in Fig. 10, a positionally adjustable photoelectric or other type of sensor 50 (having 
an operating time different from that of sensor 26A, 26B in Fig. 1) is provided at an appropriat interval L 
5 lengthwise along a belt conveyor 51 with respect to the center of primary coil 3. 

As shown in Figs. 12a and 12b, the position and operating time of sensor 50 are preadjusted so that the 
signal Q generated by the sensor is ON for the time t1 corresponding to one half (e.g. the first half) of the 
undulating waveform of filtered signals Xf and Yf generated when the object is passed through the coils by the 
conveyor 51 , beginning from the rise of the signal and continuing until the signal reaches the maximum value 
w and then falls past the zero point. 

The signals Xf and Yf are input to AND circuits 100 as shown in Fig. 1 1 together with the passage signal Q, so 
that the signals Xf and Yf can be detected only while the signal Q is ON. As shown In Figs. 12a and 12b, only the 
first halves of signals Xf and Yf are detected as signals in the form of a solitary wave similar to the weighted 
signals. 

75 Subsequently, these signals are converted to digital signals, and primary values can be obtained from the 
digital signals in the same way as for the weighted signals. 

By obtaining the primary values in this way, although there are drawbacks such as the limitation of the 
objects to packaged types and the need to adjust the position etc. of the sensor, the determination of primary 
values is simpler than by weighting. 

20 

JUDGMENT (IDENTIFICATION) — 

The MPU 15B has a program for judgement or identification as shown in Fig. 18, which will be explained later. 
Fig. 13 shows vectors representing the primary values Xp and Yp for various materials along the x and y 
axes, respectively. 

25 The phase angle of the vector for iron (Fe) is generally delayed approximately 90 degrees from that of the 
oscillation signal, and Its absolute value or output sensitivity is high. As shown in Figs. 3b and 3c, the y 
component of the signal for iron is larger In absolute value than the x component. 

In contrast, the phase difference of the vectors for stainless steel (SUS) with respect to the oscillation signal 
is not very large. For some types of stainless steel SUS1 (Fig. 13), both components are positive values. For 
30 other types of stainless steel SUS2, the x and y components are respectively positive and negative values. 
Some types of materials generate lagging signals like that of iron, and others generate leading signals. 
The material effect generally has characteristics and a phase angle which are peculiar to each Individual 
product. Although the size of the vector will change with the content volume of the product, the phase angle 
will usually not change very much. 
35 Referring to Fig. 14, points P1 (Xp1, Yp1), P2 (Xp2, Yp2), etc. represent the primary values Xp and Yp for 
various volumes of products of the same type which have material effect, but which contain no foreign matter. 

Even if products of the same type differ somewhat in volume or weight, these points will be included in an 
elliptic region or domain Dm, and will tend to be distributed near and around the central average point Po (jiX, 
jiY), which can be represented by a vector rs from the origin point 01 of the coordinate. 
40 The primary values for packages of products of the same type and substantially of the same volume and 
weight with material effect, but with no foreign matter, will fall within a small circular region Dm', with Its center 
at the average point Po, due to the vectors being substantially constant in size and angle. 

The primary values for a product having material effect and containing even minute particles of foreign 
matter, such as iron, will be represented at a point Pa (Xa f Ya) as a vector ra. Subtraction of vector rs from 
45 vector ra gives a deviation vector Ra from the average point Po. 

Thus, the region Dm defines the distribution boundary for products without foreign matter, and data such as 
point Pa outside this region are Judged to represent products which contain foreign matter. 

For products having neither material effect nor foreign matter, almost no waveform characteristics as shown 
in Figs. 3a-3c appear. In this case, random noise resulting from internal electrical factors, such as the 
50 operation of the amplifiers, wave detectors or conveyor, becomes the primary factor in determining the 
distribution region. 

As shown in Fig. 15, these products are represented as points P1, P2, etc. within a circular region Dn on the 
graph. The central average point Po of region Dn is extremely close to the origin point 01. The region Dn is 
extremely smaller than for products having material effect. 
55 In this case, the size of the deviation vector Ra is approximately the same as that of the signal vector ra. 
Because of the small size of region Dn, even If the absolute values of these vectors are relatively small, they will 
extend beyond region Dn. It Is thus possible to detect even more minute particles of foreign matter than for 
products having material effect. 
In order to obtain the data needed to determine the regions for various types of objects, there are sample 
60 tests in "M" and "N" modes, which have the following characteristics. 

M (MATERIAL EFFECT DETECTION) MODE SAMPLE TEST 

This test is carried out for objects which have material effect and contain no foreign matter by obtaining 
primary values Xp and Yp from samples as "M M mode sample data Xsm and Ysm, such as represented by 
65 points P1-P4 in Fig. 14, when the objects pass through the transducer 10. These data include the product 
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characteristics together with their occasional fluctuations, and also the effects of the various minute amounts 
of noise which remain even after filtering. 

N (NOISE DETECTION) MODE SAMPLE TEST : 

5 

This test is done for objects which have no material effect and contain no foreign matter, and for which it is 
difficult to obtain the primary values. 

When the passage of a packaged object, for example, through the transducer 1 is detected by the sensor 
26A-26B (Fig. 1), a passage signal is generated for a period of time corresponding to the time, 2t1 (Fig. 3c) 
which is the sum of time t1 preceding the passage of the center of the primary coil 3 by the object and time t1 10 
following it. While this passage signal is being generated, pairs of weighted values Xw and Yw are extracted as 
"N" mode sample data Xsn and Ysn, respectively, at intervals of a certain number of samples from among the 
continuous weighted values. For example, 10 to 20 weighted values are extracted for one packaged object. 
These values Xsn and Ysn may be represented by points P1-P4 In Fig. 15. 

Generally, because various types of random noise are a major factor concerning the signals of products 15 
which have no product characteristics, the sample values Xsn and Ysn are also random with regard to a time 
series. It is possible to appropriately extract the characteristics of the amplitude, distribution, etc. of the 
original random signals from the many data signals obtained by extracting at certain intervals as described 
above. 

The sample values Xsm and Xsn of the M and N modes are generally referred to as a sample value Xs 20 
hereinafter, and the values Ysm and Ysn as a sample value Ys. 

The sample values Xs and primary values Xp are generally termed representative values X hereinafter, and 
'Ys and Yp as representative values Y. 

CORRELATION ■ 25 

For more precise region determination, consideration is given to the correlation of representative values X 
and Y. Specifically, a correlation coefficient p as a factor representing the correlation of values X and Y is 
included in discriminant equation (5), which determines the region and will be explained later. 
The relationship between this correlation coefficient p and the region is as follows: 

Supposing X' = X - jiX and Y r « Y - *iY, in order to find a general expression, if ax and oy are selected for 30 
the unit measures of X and Y, and XVox and YVoy are used for the axes of the coordinate system in order to 
normalize it, the relationship between the coefficient p and the region can be expressed as shown in Fig. 16. 

In Fig. 16, the shape or expanse (area) of the region changes with the correlation coefficient. Specifically, if 
there is no correlation (If p - 0) between the representative values X and Y, the region is a circle, with widest 
possible width. Also, the stronger the correlation, the larger the difference between the major and minor axes, 35 
and the more elliptic the region, so that the region becomes narrower as the correlation coefficient 
approaches "1". 

if the material effect is small, because the weighted signals Xw and Yw have the characteristics of noise, that 
is, the sizes and phase of the signals are random and mutually unrelated, as shown in Fig. 15 and explained 
earlier, when the two-dimensional pairs of representative values X and Y extracted at random from weighted 40 
values Xw and Yw, respectively, are plotted on the coordinate system, they are distributed in an approximately 
circular region. 

In other words, in this case, the correlation coefficient p is a value approaching zero. 

On the other hand, if the material effect is large and the objects are either lumps or have fluctuations in their 
content volume, although the representative values X and Y which are primary values Xp and Yp will 45 
sometimes change momentarily, the ratio of them will maintain an approximately fixed relationship. In other 
words, because there is correlation between X and Y, the plotting of the pairs of these values on a coordinate 
system is generally distributed in an elliptical shape as shown in Fig. 14. 

Thus, as will be explained later, the discriminant equation which determines the region can be established as 
a discriminant equation for an elliptical region which includes correlation. By including correlation in the 50 
equation, the region can be more narrowly defined than if correlation Is not taken into consideration, and 
defective products which were mistakenly judged to be non-defective when correlation was previously not 
taken into consideration can be correctly judged to be defective, thus resulting in more accurate Judgements. 
In other words, the Introduction of the concept of correlation improves the capacity to judge the presence of 
foreign matter. 55 
. The sample test used to determine the coefficients etc. of the equation is done in the M or N mode. In this 
test, a large number of pairs of sample data Xs and Ys are obtained, stored in the memory of MPU 15B, and 
used to compute each of the following values (data for the equation). 
yX: average value of Xs 

jiY: average value of Ys 60 
ax: standard deviation regarding Xs 
oy: standard deviation regarding Ys 
oxy: covariance regarding Xs and Ys 
p: correlation coefficient expressed as 

65 
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6~X • Cy 

5 

Because the detailed explanations of the symbols listed above are contained in general reference works on 
statistics, they will be omitted here. 

As already explained, the primary values Xp and Yp are used to compute a predetermined discriminant 
equation, and judgement concerning the presence of foreign matter is made on the basis of whether or not the 
10 equation is satisfied. The equation is described below. 

Discriminant equation according to the first embodiment (basic discriminant equation) 

In the judgement processing, an equation such as the following Equation (4), which is one example of the 
statistical distribution functions found in reference works on statistics, is used as the density function of the 
15 statistical distribution of the data. This equation works both for inspection objects which have material effect 
and for those which do not, as long as the data is considered to be from the same population. 



20 f (x, y) 



27TVT :: p 2 ct v ct 



f 2 °x°y 



1 ( (X-fiX) 2 (X-MX) (Y-MY) ( Y-jttY) »' > 

: < ~ 2p + ___ L 

2( 1 7> 2 ) i V ^ ° X °y Of J 

(4) 



In this embodiment, in order to determine the region condition, the part of the above Equation (4) 
concerning the exponent of the equation is referred to as "D". This D is a numerical value which determines the 
region, and it may be determined either based on empirical data, or by the following equation which uses 
35 correlation coefficient p and the region coefficient d which will be explained later. 



D = ____ • d* 
(1 + p) 

Thus, using the above discriminant equation, the following Equation (5) is the basic discriminant equation for 
determining the region. 

1 (X-JLtX)2 (X-AtX) (Y-juY) 

1 (Y-MY) 2 2 

..... (5) 



(iy 2 ) <y i+jp 



The region coefficient "d" in the above Equation (5) is input by the region coefficient setting device 24 shown 
in Fig. 1, which is the input device for this embodiment. 

For confidence values of 2cr, 3o, etc., this coefficient d is as follows. 
60 i)For2a,d - 2, 
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2 8 
D = x d* = j ...... (6) 

1 + f> 1 + D 



ii) For3a,d = 3, 

2 18 10 
D - ; x d 2 = __ (7) 



1 + P 1 + 



This coefficient d can be appropriately selected by the operator in association with the confidence value. 
Thus, when the correlation coefficient p is calculated and the region coefficient d is selected by the operator, 
the value of "D", which determines the region, is determined. 

In addition, by using the following definitions, 

1 

A - ________ 



B = 



(1 - ji) a x « 
(1 - p 2 ) o x . 
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30 



(1 ~ ft 2 ) V 35 

Equation (5) can be rewritten as the following Equations (8) and (8)'. 
A(X-nX)* + B(X-uX)(Y-uY) + C(Y-*iY)2 « D (8) 

A(X - uX)2 + B(X -uX)(Y -uY) + C(Y - uY) 2 - D < 0 (8)' 40 

Equation (8) determines the boundary of the distribution region, and Equation (8)' is the discriminant 
equation which judges whether or not the data is included in the region. 

As explained earlier, because Equation (8) defines an elliptical or circular region, depending on the value of 
the correlation coefficient p, it is applicable to both products which have material effect and those which do 
not. 45 

Equation (8)' is included in the MPU 15B in the form of a program, and computed for each pair of primary 
values Xp and Yp at final inspection of objects. 

Theoretically, the region discriminant equation can be determined by the operator selecting and inputting 
appropriate values A, B and C, and the region coefficient d in accordance with the characteristics of the object. 
However, because It is not easy to select these values, provision has been made in this embodiment to simplify 50 
the operators selection tasks. 

Specifically, Instead of providing a device to input the values A, B and C, the MPU 15B has a program 
whereby, when the operator uses the switch 1 8 to select the equation computation mode, the data which was 
obtained and stored in the previous sample test is recalled, the statistical quantities uX, \lY, ox, ay, p, etc. are 
automatically computed, and the values A, B and C are also automatically computed. 55 

The program further computes the value D in Equation (8) with these computed results when the region 
coefficient d Is input by the operator via the device 24. This determines all coefficients of discriminant Equation 
(8)' to specify this equation. In this way, the automatic setting for the region discriminant Equation (8)' is 
completed. 

60 

Discriminant equation according to the second embodiment 

Instead of the discriminant equation for the above mbodiment, it is also possible to use the discriminant 
equation described below. The embodiment using this equation will be explained as a second embodiment 
with regard to the discriminant equation. 

if the inspection objects have an approximately constant volume or a fixed size, such as for packaged 65 
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products, even if the objects have material effect, because. the major axis of the elliptical distribution of the 
characteristics becomes shorter and the shape approaches that of a circle, the following Equations (5a) and 
(8a) can be used in place of the Equations (5) and (8). 

= D» (5a) 



0 in the first embodiment already described. 



1 



*x 2 



1 



25 Equation (5a) can be used to determine the following Equation (8a)' as the discriminant equation. 
A'(X-uX) 2 + .C'(Y-|iY) 2 -D ; < 0 (8a)' 

This case is equivalent to not taking the correlation into consideration, in other words, not using the 
correlation coefficient (p). In addition, just as in the first embodiment, it is also possible to assign a direct 
numerical value for D' using an empirical value. Discriminant equation according to the third embodiment 
30 Furthermore, instead of the discriminant equations described above for the first and second embodiments, 
it is also possible to use the discriminant equation described below. The embodiment using this discriminant 
equation will be explained as a third embodiment with regard to the discriminant equation. 

In other words, just as for the second embodiment 20 described above, if the Inspection objects are of an 
approximately constant volume or fixed size, such as packaged products, when the measurement data 
35 (primary values) P1-P4, etc. of a typical inspection object are plotted on an x-y coordinate system, they are 
gathered together in the vicinity of a specific location on the coordinate system, as shown in Fig. 17. 

These data can be considered to be distributed within a closed range of predetermined spreads of ±d.ox 
(= ± Dx) and ±d. ay(- ±Dy) in the positive and negative directions with respect to the average values uX 
and uY along the x and y axes, respectively. This derives the following Equations (5b) and (5c). 
40 uX-Dx ^ X < uX + D x (5b) 
. uY-D y < Y ^ uY + D y (5c) 

The above Equation (5b) and (5c) can be transformed to obtain the following Equations (8b)' and (8c)'. 
IX-uX|-Dx £ 0 (8b)' 
. |Y-uY|-Dy ^ 0 (8c)' 
45 Equations (8b)' and (8c)' are the discriminant equations, and objects which satisfy both of these equations 
are non-defective products containing no foreign matter. On the other hand, objects which fail to satisfy either 
one or both of these equations are judged to be defective products containing foreign matter. An illustration of 
the region determined by these equations appears as a region enclosed within a rectangle, as shown in Fig. 17, 
In other words, non-defective and defective products fall respectively in and outside the region. 
50 In Fig. 17, the distance between the central point Po of the region and the origin point 01 of the coordinate 
system represents the average size of the material effect of the product. Thus, for objects having no material 
effect, these points Po and 01 are the same. 

Discriminant equations (8b)' and (8c)' are also programmed in the control device 15. The correlation 
coefficient p is also not taken into consideration for this third embodiment. The values Dx and Dy may also be 
55 empirical values. . 

As explained above, the operator operates the following devices: . 
the device 24 to specify the region coefficient d; 

. the numeric keys of keyboard 25 (Fig. 1) to specify the code number for the type of object; 
the command key Mi of keyboard 25 to input and store in the MPU 15B with the code number the various 

60. coefficients, such as A-D for the first embodiment, A', C and D' for the second embodiment, or Dx and Dy for 
the third embodiment, and the statistical quantities uX and uY of the discriminant equation/s computed with 
the data obtained from the sample test results and with the coefficient d input by the operator; and 
the command key Mo to r call with the code number the data when the same discriminant equation is to be 
used on a subsequent occasion. 

65 Through these operations, various coefficients, statistical quantities, tc. related to each of a plurality of 
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This Equation (5a) is applicable when p = 
Here, with the following definitions, 



15 A' = 



20 
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different product types are stored during the first inspection, and, if a product of one of these types is to be 
inspected again, the relevant information can be recalled and instantly used. 

OPERATION 

The following is an explanation of the procedures involved in the detection of foreign matter. 5 
PRELIMINARY PROCEDURES 
(1) Test Operation 

The operator sets the selector switch 18 to the sample test mode. At the same time, the RAp indicator 20 is 10 
set to indicate the absolute value r and phase angle 8 of a vector r from the origin point on the co ordinate 
system (Fig. 14, 15). The values r and 0 are calculated with weighted values Xw and Yw as follows, r « j/Xw 2 + 
Yw* 



7L XW 



75 



In this way, the operator is able to continuously monitor the Indication values during no-load operation and 20 
the subsequent sample test. 

First, a no-load test operation is performed to find the noise level of the apparatus Itself on the indicator 20. 

"Next, the operator selects a number of samples, which have been determined to be freeof foreign matter, 
from a type of objects to be inspected. He puts the samples to be passed through the coils, reads the absolute 
values on the indicator 20, and determines whetheror not the values have increased from those during the 25 
no-load operation. 

If the values are clearly larger than the noise level, it is judged that there is material effect, and the 
subsequent sample test is carried out in the M mode. If not, it is judged that there is no materia! effect, and the 
sample test is carried out In the N mode. 

30 

(2) Sample Test 

The operator sets the selector switch 23 to either the M or N mode, presses the execution button 21 , selects 
a large number of samples having no foreign matter from that type of objects, and puts the samples to be 
passed through the coils. As a result, regardless of which mode was used, the data for the samples is stored in 
the MPU 15B. 35 

(3) Discriminant Equation Determination (Specifying) 

When the operator sets the switch 18 to the equation computation mode and presses the execution button 
21 , the values pX, jiY, ax, ay, p, etc. for the samples are computed as described earlier with the data obtained 
in the sample test. - 40 

For the first embodiment, the equations for the values A, B and C are then calculated to determine the 
coefficients for the left side of Equation (8). 

When the operator inputs the appropriate region coefficient d, the value D which determines the region 
conditions is calculated, and the discriminant equation (8)' is specified. This completes the preparations for 
the computation of the discriminant equation for the first embodiment. 45 

For the second and third embodiments as well, the equations (8a)' or (8b)' and (8c)' are determined 
likewise. 

The values D, D' Dx, Dy, etc. can be found by calculating as explained, or, if empirical values are known, they 
can be input directly. 

50 

(4) Storage of Coefficients and Statistical Quantities 

When the operator completes the equation determination, it is possible, Instead of proceeding directly to 
actual operation, to store the coefficients and statistical quantities for the discriminant equation peculiar to 
each type of objects and recall them at any time, as explained earlier. 

Specifically, the operator stores the values A, B. C, D, pX, jiY, etc. for the particular discriminant equation by 55 
inputting the code number for the type of objects via the keyboard 25 and then pressing its memory key Mi, 
and recalls the data when needed by Inputting the same code number and then pressing the key Mo. 

Thus, in actual operation, by specifying a code number for various data pertaining to a particular type of 
object and storing it, part of the preliminary procedures can be eliminated during subsequent inspection work. 

60 

ACTUAL OPERATION 

Actual operation begins once the preliminary procedur s have been completed. During actual operation, 
the R/<p indicator 20 is set to indicate the absolut value R and phas angle <p of a deviation vector R, which are 
computed by Equations (11) and (12) below with the primary values Xp and Yp. 

When the operator sets the switch 1 8 to the operation mode for actual operation, the program control circuit 65 
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17 inputs the program for inspection of objects to the MPU 15B. In this mode, the judgement section of MPU 
15B receives the data in the M mode. 

The operator presses the execution button 21 , and passes through the transducer 10 objects the samples 
of which have been tested, but it is unknown whether or not the objects contain foreign matter.. If the 
5 transducer 10 outputs a signal XY as shown in Fig. 3a for one object, it is processed through th various 
devices and input to the MPU 15B. 

In the MPU 15B, a pair of primary values Xp and Yp are obtained by the program of Fig. 8. Even during this 
processing, other objects are continuously detected (parallel processing). . 

As shown In Fig. 18, the following equa tions are then computed w ith the primary values, and the results are 
10 indicated in the RAp indicator 20. R _ j/(Xp = u*) 2 + (Yp-pY) 2 (11) 



15 



35 



As mentioned earlier, the values uX and jiY are the average sample values obtained In the sample test in the 
M or N mode, and the values A-D, or A', C and D', or Dx and Dy in the discriminant equation/s (8)', (8a)' f or 
20 (8b)' and (8c)' have been predetermined in the preliminary procedures. 

The equation/s Is/are then computed with the values Xp and Yp substituted for X and Y, respectively. If the 
result satisfied the equation (i.e., the left side is 0 or negative), the data is judged to be in the region. If the 
equation is not satisfied, the data is judged to be outside of the region. . 

If the data is outside the region, the alarm 19 is activated to warn that the object has not passed the 
25 inspection, and a sorting device (not shown) further down the line is commanded to remove the object. 

With this method, when an object contains water, salt or the like and has a material effect, if the amount of 
the contained material varies, the electromagnetic characteristics will change and the region conditions of the 
reference object (standard product) will not be satisfied. This will cause a signal which indicates that the object 
is defective to be generated. Thus, this method may be effective for use in quality control. 
30 The RAp indicator 20 may be adapted to indicate the type of material effect of the object, or the type of 
foreign matter (metal, iron or non-ferrous metal, etc.) which has been detected in the object, as shown in 
Fig. 13. 



ADVANTAGES OF DETECTOR WITH JUDGEMENT SECTION 



(1) Because the data obtained in the sample test is automatically computed to determine the 
judgement conditions, the determination of the judgement conditions is almost completely automatic. In 
addition, because the data is stored after it is collected and can then be recalled and used again for the 
, same inspection object, there is no need to collect the data for the inspection object at that time. 
40 (2) The operator only has to consider such things as whether to select a confidence level of 2a or 3a In 

order to determine the region coefficient, and there is no need for adjustments which require skilled 
expertise, such as phase adjustments of the excitation signal or signal detection reference signal, as are 
required in prior art detectors. 

(3) It Is possible to achieve a remarkable improvement over prior art detectors in the precision of the 
45 detection of foreign matter which is contained in Inspection objects having material effects. 

(4) It is possible to inspect more types of inspection objects than with prior art foreign matter detectors 
of this type, thus Increasing the range of Inspection applicability. In other words, regardless of whether 
the inspection objects are separate or packaged, and also regardless of whether they have material effect 
or not, the detector of this Invention can be set and used according to a fixed procedure without being 

50 affected by the configuration or characteristics of the inspection objects. 

General Advantages of the Invention 
With this invention, as described above, it is possible to perform accurate detection even if the inspection 

object has material effect, and also even if the inspection object contains minute particules of metal which 
55 react in the same phase as that material effect. 

Moreover, by knowing in advance the relationship between the type of metal and the phase, it is also 

possible to detect the type and size (quantity) of the metal. Also, because there Is no need for precise 

trial-and-error adjustments of either the phase or the sensitivity, as was required in prior art detectors, this . 

detector can be easily operated by anyone without requiring extended preparations. 
60 Thus, damage to processing equipment which might result from the presence of foreign matter in raw 

materials, etc. is prevented, and, furthermor , the shipment of products etc., containing foreign matter is also 

prevented. 

In addition, because the detection signals contain not only data on whether or not foreign matter Is present, 
but also quantitative data as well, the detection of foreign matter or of quality differences is easily 
65 accomplished. Thus, with this detector, it is also possible to detect instances where the quality of the product 
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differs from the specified level. 

Furthermore, because the composition is simple, the detector can be provided at a lower cost than the prior 
detectors described herein, and maintenance and adjustments can be easily performed. 

Even if metal or other foreign matter is contained inside an inspection object having material effect and that 
foreign matter has a phase of the same direction as the phase direction of said material effect, because the 5 
two-dimensional pair of digital primary value signals for the inspection object is computed and identified using 
a discriminant equation which has been predetermined in consideration of that material effect of the inspection 
object, identification can be performed with greater sensitivity than possible using a detector of the prior art. 

For the same reason as above, even if there is noise generated by the device, the detection of metals and 
other foreign matter can be performed without being affected by that noise. 10 

Because the detection signals are expressed by two coordinates (i.e., in a two-dimensional plane) the type 
and size of metal can be easily identified, and even foreign matter which generates no more than a small signal 
can be detected. 

15 

Claims 

1 . A detector comprising: 

means (1 to 10) for providing a two-dimensional pair of analog signals representative of electromagnetic 
parameters of an object; 20 
means (11,12) for converting said analog signals into respectively a first and a second series of digital 
values; 

' means (14,15) responsive tc said "two series of digital vaicies for determining respective representative 
values X and Y; and 

computation means (15) for performing a computation based on the following discriminant equation(s) (i), 25 
or (II) in order to determine whether the conditions expressed thereby are met: 

(1) A' (X - uX) 2 + C'(Y - u.Y) 2 - D' < 0 
where the coefficients A', C, D', uX, and u,Y are values approximately selected for the object; 

(H) |X-uX| -Dx<0and 

IY-uYI - Dy £ 0 30 
where the coefficients uJ(, uY, Dx and Dy are values appropriately selected for the object. 

2. A detector as claimed in claim 1, wherein the converting means (11,12) is capable of providing a 
plurality of sets of first and second digital values, each set being produced in response to analog signals 
representative of electromagnetic parameters of a respective sample of said object; wherein, with regard 

to equation (i): 35 

1 

A' - . . 

ax 2 

40 

1 •• 

. C' = _____ : 

ay 2 

/. 45 

D' « 2d 2 ; wherein, with regard to equation (H): 
Dx — d.ox 
Dy « d.oy 
and wherein: 

uX and u,Y are averages of representative values X and Y, respectively, for the samples ; and 50 
ox and ay are respectively standard deviation values of X and Y for the samples; 

the detector further comprising means (15) for obtaining the values of said coefficients uX and jiY, 
together with A', C and D' or Dx and Dy, from said sets of first and second digital values and a value d 
appropriately selected for said object. 

3. A detector as claimed in claim 2, including means (15,25) permitting storage of data defining 55 
coefficients to be used in a discriminant equation together with a code associated with the type of object 

to which the coefficients relate, thereby permitting later retrieval of said data by entry of said code 
number. 

4. A detector comprising: 

means (1 to 10) for providing a two-dimensional pair of analog signals representative of electromagnetic 60 
parameters of an object; 

means (11,12) for conv rting said analog signals into respectively a first and a second series of digital 
values; 

means (14,15,50) for deriving from said first and second series of digital values respective representative 
values X and Y, said deriving means including means (50) for generating a passage signal while said 65 
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object is passing through a predetermined region during which said analog signals are provided, and 
means (15) for selecting said representative values X and Y from said two series of digital values 
produced while said passage signal is generated; and 

comparing means (15) for comparing said representative values X and Y selected while said passage 
5 signal is generated with predetermined reference values. 

5. A detector as claimed in claim 4, wherein said providing means (1 to 10) is so arranged that the 
analog signal each adopt values which differ in a first sense from an initial value as the object travels 
through a first part of its path and which differ in a second, opposite sense from the initial value as the 
object travels through a second part of its path, and wherein said predetermined region corresponds to 

10 one of said parts of the object's path. 

6. A detector as claimed in claim 4 or 5, including means (15,25) for storing a plurality of sets of 
reference values, each set relating to a respective type of object, said means (15,25) permitting retrieval 
of each set in response to entry of a respective code. 

7. A detector comprising: 

15 means (1 to 10) for providing a two-dimensional pair of analog signals representative of electromagnetic 

parameters of an object; 

means (11,12) for converting said analog signals into respectively a first and a second series of digital 
values; 

means (14,15) for deriving from said first and second series of digital values respective representative 
20 valuesXand Y; 

• means (15,25) permitting storage and retrieval of a selected cne-of n plurality of.stored sets .of reference 
values, each set being associated with a respective type of object; and 

comparing means for comparing the representative values X and Y with the retrieved set of reference 
values. 

25 8. A detector according to any one of claims 4, to 7, wherein said comparing means (15) is operable to 

compare the representative values X and Y with reference values by determining whether the conditions 
expressed by one or more discriminant equation(s) are met. 

9. A detector according to claim 8, wherein said equation is: 
A(X - pX) 2 + B(X - |iX)(Y - pY) + C(Y - pY) 2 - D ^ 0 

30 where the coefficients A, B, C, D, pX and jiY are values appropriately selected for the object. 

10. A detector according to claim 8 t wherein said equation is: 
A'(X - pX)2 + C'(Y - p.Y)2 - D' £ 0 

where the coefficients A f , C, D', pX, and pY are values appropriately selected for the object 

1 1 . A detector according to claim 8, wherein said equations are : 
35 |X-pX|-Dx<0 

IY-|iY|-Dy>u 

where the coefficients pX, pY, Dx and Dy are values appropriately selected for the object. 

12. A detector comprising: ■ 

means (1 to 10) for providing a two-dimensional pair of analog signals representative of electromagnetic 
40 parameters of an object; 

means (11,12) for converting said analog signals into respectively a first and a second series of digital 
values; 

means (14,15) responsive to said two series of digital values for determining respective representative 
values X and Y; and 

45 computation means (15) for performing a computation based on the following discriminant equation(s) (i), 

(ii) or (iii) in order to determine whether the conditions expressed thereby are met: 

(i) A(x - pc)2 + B(X - pX)(Y - p.Y) + C(Y - jiY)2 - D ^ 0; 

(ii) A'(X - pX) 2 + C'(Y - pY) 2 - D' £ 0; 

(iii) |X-pX| -Dx^Oand 
50 |Y- jxY| -Dy^O; 

the detector futher comprising means (15,25) permitting storage of the coefficients A, B, C, D, pX, pY, or 
A' f C\ D', pX and jiY, or pX, pY, Dx and Dy and subsequent retrieval thereof for use by said computation 
means (15). 

13. A detector as claimed in claim 12, wherein the converting means (11,12) is capable of providing a 
55 plurality of sets of first and second digital values, each set being produced in response to analog signals 

representative of electromagnetic parameters of a respective sample of said object; . 
wherein, with regard to equation (i): 
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( l-jp 2 )cx* . 
)axay 

l /o 



(1-J> 2 )ay 
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D = ' .d 2 , 

(1+JP) 

wherein, with regard to equation (ii) : 
.A. I = 



is 



20 

1 



(XX* 

25 



ay 2 

D' = 2d 2 ; 



5(7 



wherein, with regard to equation (iii) : 
Dx = d.ox 

Dy = d.ay 35 
and wherein: 

joX and \lY are averages of representative values X and Y, respectively, for the samples; 
ax and ay are respectively standard deviation values of X and Yfor the samples; and 
p is a correlation coefficient of X and Y for the samples ; 

the detector further comprising means (15) for obtaining the values of said coefficients jiX and ^lY, 40 
together with coefficients A, B, C, and D, or A', C and D', or Dx and Dy. from said sets of first and second 
digital values and a value d appropriately selected for said ofaject, thereby permitting subsequent storage 
of said coefficients. 

14. A detector as claimed in claim 12 or 13, wherein said means (15,25) permitting storage is arranged to 
allow a plurality of sets of coefficients, each set being associated with a respective type of object, to be 45 
stored and subsequently selectively retrieved. 

1 5. A detector according to any preceding claim, wherein the providing means (1 to 10) comprises: 
an oscillator (1) for generating a basic oscillating signal; 

a phase shifter (8) connected to said oscillator (1) and adapted to provide a phase-shifted oscillating 
signal; 50 
an electromagnetic transducer (10) including an excitation coil (3) connected to said oscillator (1), and 
two interconnected detection coils (4,5) magnetically coupled to said excitation coil and adapted to 
produce a differential signal therebetween when an object is passed through said transducer (1Q) ; 
a first detector (7A) connected to said detection coils (4,5) and said oscillator (1 ) and adapted to produce 
a first detected analog signal representing the component of said differential signal in phase with said 55 
basic oscillating signal; and ; 

a second detector (7B) connected to said detection coil (4,5) and said phase shifter (8) and adapted to 
produce a second detected analog signal representing the component of said differential signal in phase 
with said phase-shifted oscillating signal. 

.' • •' 60 
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